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ABSTRACT 



O: 
U: 

O ', Clusters of galaxies are uniquely important cosmological probes of the evolution 

f-) . of the large scale structure, whose diagnostic power depends quite significantly on the 

ability to reliably determine their masses. Clusters are typically modeled as spherical 
systems whose intracluster gas is in strict hydrostatic equilibrium (i.e., the equilib- 
rium gas pressure is provided entirely by thermal pressure), with the gravitational 
psj ' field dominated by dark matter, assumptions that are only rough approximations. In 

, fact, numerical simulations indicate that galaxy clusters are typically triaxial, rather 

■ than spherical, and that turbulent gas motions (induced during hierarchical merger 

events) provide an appreciable pressure component. Extending our previous work, we 

OA , present results of a joint analysis of X-ray, weak and strong lensing measurements of 

Abell 1689. The quality of the data allows us to determine both the triaxial shape of 
£f) , the cluster and the level of non-thermal pressure that is required if the intracluster gas 

■ is in hydrostatic equilibrium. We find that the dark matter axis ratios are 1.24 ± 0.13 
' and 2.02 ± 0.01 on the plane of the sky and along the line of sight, respectively, and 

that about 20% of the pressure is non-thermal. Our treatment demonstrates that 

■ the dynamical properties of clusters can be determined in a (mostly) bias-free way, 
• rH ' enhancing the use of clusters as more precise cosmological probes. 

X ■ 

5^ , Key words: cosmology: observations - galaxies: clusters: general - galaxies: clusters: 

individual (Abell 1689) - gravitational lensing: strong - gravitational lensing: weak - 
X-rays: galaxies: clusters 



1 INTRODUCTION 

Clusters of galaxies are the largest bound systems, formed at 
relatively late times. As such, their mass function sensitively 
depends on the evolution of the large scale structure (LSS) 
and on the basic cosmological parameters. Given this great 
potential, the use of clusters as a cosmological probe hinges 
on our ability to accurately determine their masses. 

Mass determinations based on X-ray observations cus- 
tomarily assume spherical symmetry and strict hydrostatic 
equilibrium (HE), i.e., that intracluster (IC) gas pressure is 
provided entirely by thermal motions, P to t = -Pth- Under 
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these assumptions X-ray measure ments can be s uccessfully 
used to constrain the mass profile (|Sarazinl ll988h 

In addition to this X-ray based method, optical obser- 
vations of giant arcs and multiple images produced by strong 
gravitational lensing (SL) in the central parts of clusters and 
slight distortions of background sources in the weak lensing 
(WL) regime allow us to determine projected cluster masses 
without invoking the assumptio n of hydrostatic equilibrium 
(|Miralda-Escude fe Babull995n . but an a priori assumption 
about the spherical symmetry is still often used to deduce 
the three-dimensional mass profiles. 

However, N-body simulations indicate that dark mat- 
ter (DM) halos are triaxial with intermediate-major and 
minor-in termediate axis ratio s typically of the order of 
~ 0.8 (|Wang fe White! l2009h . and hydrodynamical nu- 
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merical simulations suggest that even after equilibrium 
is established a significant fraction of the pressure sup- 
port agai nst gravity comes from subsonic non-thermal gas 
motions dLau et al.ll2009l; IZhang et al.l [20101 ; iRichard et al.l 
l2010l : lMeneghetti et al.ll2010l ). raising doubts on the viability 
of cluster sphericity and purely thermal gas pressure. 

Moreover, since lensing is sensitive to the integrated 
mass contrast along the li ne of sight, both dep artures from 
the spherical assumption ()Morandi et a .1 201dh and the in - 
clusion of non-thermal pressure support ( Molnar et al.ll20icil ) 
can explain the long-standing discrepancy between cluster 
masses determined from X-ray and gravitational lensing ob- 
servations, the latter often being significantly higher than 
the former. 

The main goal of our work is to resolve these dis- 
crepancies between X-ray an d gravitational lensing mass 
for the the cluster Abell 1689 dAndersson fc Madeiskill2004l; 
Lemze et a l. 2008; Ricrncr-S0 rensen et al .1120091 ; iPeng et al.l 
20091 : iMorandi et al.ll201lh . We tackle this by developing a 



novel method in order to infer both the desired triaxial shape 
and physical properties of clusters and the non-thermal pres- 
sure via a joint X-ray, weak and strong lensing analysis. 

We ex tend the findings of our previous work on 
Abell 1689 |Morandi et al.ll201ll ), where we combined X-ray 
and strong lensing data; in that work we accounted for the 
three-dimensional geometry, which allowed us to resolve the 
discrepancy between the mass determined from X-ray and 
strong gravitational lensing observations in the inner region 
(R <J 400 kpc), assuming strict HE. In the present paper 
we jointly analyze also weak lensing data that map the pro- 
jected mass profile out to ~ 3 Mpc. With this additional 
data we are also able to determine the non-thermal pressure 
of the IC gas. 

Hereafter we assume the flat ACDM model, with mat- 
ter density parameter Q m = 0.3, cosmological constant 
density parameter Qa = 0.7, and Hubble constant Ho = 
70 km s _1 Mpc -1 . Unless otherwise stated, quoted errors are 
at the 68.3% confidence level. 



2 DATASETS AND ANALYSIS 

A full description of the X-ray and SL analysis can be found 
in IMorandi et all (12011); for details on the WL analysis we 
refer to Umetsu et al. (2009). Here we only briefly summa- 
rize the most relevant or novel aspects of our data reduction 
and analysis of Abell 1689. 



2.1 X-ray analysis 

We analyzed two sets of Chandra observations (ID numbers 
6930 and 7289) from the NASA HEASARC archive with a 
total exposure t ime o f approximately 150 ks. With respect to 
IMorandi et al.l |201ll ) , we repeated the X-ray data reduction 
and analysis by implementing the most recent Chandra cal- 
ibrations: we used the CIAO software (version 4.3) and the 
gain file provided within CALDB (version 4.4.3). We mea- 
sured the gas density profile in a non-parametric way from 
the surface brightness recovered by a spatial analysis, and 
we inferred the projected temperature profile by analyzing 
the spectral data. 

The X-ray images were extracted from the two event 



files in the energy range 0.5 — 5.0 keV, then corrected by the 
exposure map to remove the vignetting effects, followed by 
the masking out of point sources. We constructed a set of 
n = 57 elliptical annuli of minor radius r m around the cen- 
troid of the surface brightness with eccentricity e b i (r) fixed 
to that predicted from the eccentricity ey (r Jj of the DM 
halo on the plane of the sky as determ i ned fr om the SL 
data (see Section [2~4l and IMorandi et all (|2010l )). We then 
deduced the electron density n e — n e (r;e c i) by deproject- 
ing the surface brightness profile, obtaining n = 57 radial 
measurements in ellipsoidal shells. Note the dependence of 
n e (r; e c /) on the eccentricity e c i of the IC gas along the line 
of sight, still to be deter mined (for further de tails see below 
and also Appendix A of lMorandi et all boiO). 

The spectral analysis was performed by extracting the 
source spectra from n* (n* = 9) elliptical annuli of minor 
radius r* m around the centroid of the surface brightness and 
with eccentricity as predicted from that of the the DM halo 
(as above). For each of the n* annuli the spectrum was ana- 
lyzed by simultaneously fitting absorbed MEKAL models to 
the two observations, in the energy range 0.6-7 keV (0.9-7 
keV for the outermost annulus only) ; we fixed the redshift to 
z = 0.183 and the photoelectric absorption to the Galactic 
value. For each of the annuli, we considered three free pa- 
rameters in the spectral analysis: the normalization of the 
thermal spectrum K\ oc J dV, the emission-weighted tem- 
perature Tproj i, and the metallicity Z\. 



2.2 Strong lensing analysis 

For the strong l ensing analysis we refer to the findings of 
iLimousin et al.l (|2007l ). who presented a reconstruction of 
the mass distribution of the galaxy cluster Abell 1689 using 
detected strong lensing features from deep ACS observations 
and extensive ground based spectroscopy. They presented a 
parametric strong lensing mass reconstruction using 34 mul- 
tiply imaged systems, and they inferred two large-scale dark 
matter clumps, one associated with the center of the cluster 
and the other with a northeastern substructure. We masked 
out the north-eastern sector of both the 2D projected mass 
map and the X-ray data, in order to avoid the contribu- 
tion from this secondary substructure. We masked out also 
the central 25 kpc, which is affected by the mass distribu- 
tion of the cD galaxy. From the strong lensing analysis the 
major clump of the cluster appears to be elongated with a 
major-minor axial ratio on the plane of the sky of 1.24±0.13 
and a position angle of 0.4 ± 1 degrees. We rebinned the 
two-dimensional projected mass map into elliptical annuli, 
whose eccentricity, cen troid and position ang le are the same 
as those inferred from ILimousin et all l|2007l ). Then we cal- 
culated average values of the elliptical projected mass profile 
E(_R), R being the minor radius of the two-dimensional el- 
liptical annuli. We also calculated the covariance matrix C s i 
among all the measurements of £(i?) similarly to the WL 
analysis (see H2.3l for further details). 



1 ey = \/T— (b' /a') 2 , a'(b') being the major (minor) axis on the 
plane of the sky. 
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2.3 Weak lensing analysis 



2.4 Joint X-ray+Lensing analysis 



Fo r the weak le n sing analysis we refer to the findings 
of lUmetsu et alJ ([2009) , who derived a projected two- 
dimensional mass map of Abell 1689 based on an entropy- 
regularized maximum likelihood combination of lens magni- 
fication with distortion of red background galaxies in deep 
Subaru images. Advantages of that approach are: 1) by com- 
bining the distortion and magnification measurements the 
convergence can be obtained unambiguously with the cor- 
rect mass normalization, while the convergence derived from 
distortion data alone suffers from the mass-sheet degener- 
acy; 2) the method is not restricted to the weak regime 
but applies to the whole area outside the tangential critical 
curve, where nonlinearity between the surface mass density 
and the observables extends out to a radius of a few arcmin- 
utes. In the two-dimensional projected mass map in WL we 
masked out the region corresponding to the secondary sub- 
structure seen in the SL regime ( i|2.2) l. and then rebinned 
the data into elliptical annuli, whose eccentricity, centroid, 
and position angle are the same as that from SL. 

In order to calculate average values of the elliptically- 
symmetric projected mass profile E(P), where R is the 
minor radius of the two-dimensional elliptical annuli, we 
compute E = E(P) from a weighted radial projection 
of the two-dimensional sur face mass density S. Following 
iHobson fc Maisinged l|2002ft , we can formally express a lin- 
ear relation between £ and E: 



-M(r) E + Af, 



(1) 



where TV is the noise vector (at the i'th pixel) with covari- 
ance matrix C, and M = M(W) is the lD-to-2D mapping 
matrix. M is defined as an A r p i x x iV ann matrix, with N p i* the 
total number of pixels, iV ann the total number of elliptical an- 
nuli, and elements M.™ = Wi m , where Wi m (0 ^ Wi m < 1) 
is the fraction of the area of the i'th pixel lying within the 
m'th elliptical bin. An optimal solution for this problem is 
given by: 



E = {M'C^My^'C 1 s 



(2) 



The binned surface mass data E have covariance matrix 



Cwi = {M'C^MY 



(3) 



We use Monte Carlo integration to calculate Wim, while we 
refer E to th e area-weighted minor radius R m of the m'th 

elliptical bin (|Umetsu fc Broadhurst 2008) : 

Umctsu & Broadhurst] IMPOST ), who 



Note that unlike 



used 2D data in order to retrieve best-fit model parame- 
ters, in the present paper we used rebinned surface mass 
profiles. While we verified that this does not change signif- 
icantly the desired best-fit model parameter ( t|2.4p . fitting 
rebinned data has the advantage of requiring significant less 
computational time and it is usually a better approach when 
it comes to deal with real data, such as WL, possibly affected 
by systematics. Indeed rebinning involves a smearing of the 
data and also of their errors (statistical + systematics) on a 
larger area: provided that the systematics can be regarded 
as symmetric errors on a large enough area, this is usually 
a more robust approach than relying on 2D data. 



Here we bri e fly s ummarize the major findings of 
iMorandi et al.l (|2010l ) for the joint X-ray+Lensing analysis 
in order to infer triax ial physical propertie s : add itional de- 
tails can be found in IMorandi et all (|2007l , l20l6T ). We also 
outline the improvements implemented in the present work, 
i.e., the inclusion of WL data and the measurement of the 
non-thermal component of the IC gas. 

The lensing and the X-ray emission both depend on the 
properties of the DM gravitational potential well, the for- 
mer being a direct probe of the two-dimensional mass pro- 
file and the latter an indirect proxy of the three-dimensional 
mass profile through the hydrostatic equilibrium equation 
applied to the gas temperature and density. In order to in- 
fer the model parameters of both the IC gas and of the 
underlying DM density profile, we perform a joint analysis 
of SL, WL and X-ray data. We briefly outline the methodol- 
ogy for inferring physical properties in triaxial galaxy clus- 
ters: (1) We start with a generalized Navarro, Frenk and 
White (gNFW) tri axial model of the DM as described in 
IJing fc Sutd (|2002l ), which represents the total underlying 
mass distribution and depends on a few parameters to be 
determined, namely the concentration parameter c, the scale 
radius r s , the inne r slope of the DM a and the two axis ra- 
tios; (2) following ITee fc Sutd (120031 . 120041 ). we recover the 
gravitational potential and surface mass profile E of a dark 
halo with such a triaxial density profile; (3) we solve the gen- 
eralized hydrostatic equilibrium equation, i.e., including the 
non-thermal pressure (Equation Q below), in order to infer 
the theoretical three-dimensional temperature profile T gas in 
a non-parametric way, given IC gas with a density profile as 
deprojected from the X-ray data (see §2.1) and sitting in the 
gravitational potential well previously calculated; finally, (4) 
the joint comparison of T gas with the observed temperature 
and of E with the observed surface mass density gives us 
the parameters of the triaxial DM density model and the 
non-thermal component of the gas, and therefore all the de- 
sired physical properties of the IC gas and the DM triaxial 
ellipsoids. 

In particular, for the X-ray analysis we rely on a gener- 
alizat i on of the hydrostati c equilibrium equation (|Lau et al.l 
l2009l : iMolnar et all |2010h . which accounts for the non- 
thermal pressure P n t and reads: 

VPtot = -p s ^V(t> , (4) 

where p gas is the gas mass density, (j> is the gravitational 
potential, and P to t = Pth + -Pnt • 

In this first implementation of our new approach, and 
given the limited ability to carry out a more detailed spatial 
analysis, we model P n t by assuming that the non-thermal 
pressure of the gas is a constant fraction £ of the total pres- 
sure Ptot, i.e., 



(5) 



From Equations @ and ([SJ we point out that neglecting 
Pnt (i.e., setting P to t = Pth) systematically biases low the 
determination of cluster masses, roughly by a factor of £. 
Note that X-ray data probe only the thermal component of 
the gas P th = n e kT gi 



The work of iLee fc Sutcl l|2003h showed that the IC gas 
and DM halos are well approximated by a sequence of con- 
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centric triaxial distributions with different eccentricity ratio. 
We define ey (ey) and ey (e c i) as the eccentricity of DM 
(IC gas) on the plane of the sky and along the line of sight, 
respectively. The iso-potential surfaces of the triaxial dark 
halo coincide also with the iso-density (pressure, tempera- 
ture) surfaces of the intracluster gas. Th is is simply a direct 
conse quence of the X-ray shape theorem (|Buote fc Canizaresl 
19!) 1 :•: the hydrostatic equilibrium equation Q yields 

VP x X7<j) = Vp x V0 = 0. (6) 

Note that ey — ey (ey , u, a) and ey = ey (ey, u, a), 
with u = r/r s , unlike the constant ey,ey for the adopted 
DM halo profile. In the whole range of u, ey /ey (ey /ey) is 
less than unity (~ 0.7 at the center), i.e., the intracluster gas 
is altogether more spherical than the underlying DM halo. 

We construct the likelihood performing a joint analysis 
for lensing and X-ray data, in order to constrain the prop- 
erties of the model parameters q 

q= (c,r s ,a,e c /,£) (7) 

representing the concentration parameter, scale radius, in- 
ner slope of the DM, eccentricity of the DM along the line 
of sight and fractional non-thermal pressure of the gas, re- 
spectively. 

Therefore, the system of equation we simultaneously 
rely on in our joint X-ray+Lensing analysis is: 

Tgas = Tg as (c,r s ,a,e c ,,£) (8a) 

E = E(c,r s ,a,e c ,) (8b) 

where the three-dimensional model temperature T gas is re- 
covered by solving equation Q and constrained by the 
observed temperature profile, and the model surface mass 
profile E is recovered by projection of the triaxial gNFW 
DM model and constrained by weak/strong lensing mea- 
surements. 

The method works by constructing a joint X- 
ray+Lensing likelihood: 

£ — £x ■ £lcns = £x ■ £sl ■ £wl , (9) 

£x and £i cn s being the likelihoods coming from the X- 
ray and lensing (WL+SL) data, respectively. Now, £x oc 
exp(-Xx/2), where 

Xx = E 2 (10) 

i-1 a T* ■ 

z — 1 proj,i 

Tproj.i being the observed projected temperature profile 
in the z't h ring and r pr oi.i(q) the convenient projection 
(following Ma zzotta et al.l 12004 ) of the theoretical three- 
dimensional temperature T gas (q); the latter is the result 
of solving the hydrostatic equilibrium equation, with the 
gas density n e (r; ey) inferred from the X-ray surface bright- 
ness, and assuming a gNFW parametrization for the DM 
POM = pDM{r;c,r s ,a,ey). The lensing contribution is 

_ exp{-i[S-S'] t C7- 1 [S-S1} nu 

lmS (27T) m */2|C| 1 /2 ' ^ 

where C is the covariance matrix of the projected mass pro- 
file from strong (C s i) and weak (C w i) lensing data, |C| in- 
dicates the determinant of C, E* = (E*, E|, E^«) are 
the observed measurements of the projected mass profile in 



the m* elliptical annuli, and E(c, r s , a, ey) is the theoretical 
projected mass profile within our triaxial DM model. 

For the covariance matrix of the strong lensing measure- 
ments, we also account for systematic errors by replacing the 
observed covariance matrix with the following expression: 

Csl > Csl + CTgys X (12) 

where I is the identity matrix. Since we are using the sur- 
fa ce mass density as outp ut from the parametric SL analysis 
imousin et all (|2007D , we introduce the parameter o~ S y S 
to allow some extra freedom due to possible systematic bias 
resulting from their analysis. In order to calculate <7 ay s we 
assumed that systematic errors can be described as Gaus- 
sian errors via a diagonal matrix a"g ya ^ with the same value 
in each of the diagonal elements. We checked that this sim- 
plified assumption does not significantly affect the average 
value of the physical parameters, while it slightly increases 
(12%-20%) their errors with respect to the case where we 
neglect ovys. 

We evaluated the probability distribution function of 
model parameters £ = £(q, ovys) via the Markov Chain 
Monte Carlo (MCMC) algorithm. Thus we can determine 
our best-fit model parameters and their errors, and derive 
various properties of the cluster. 



3 RESULTS AND DISCUSSION 

In the previous section we showed how we can determine 
the physical parameters of the cluster by fitting the available 
data based on the hydrostatic equilibrium equation and on a 
DM model that is based on robust results of hydrodynamical 
cluster simulations. In this section we present our results and 
discuss their main implications. We particularly focus on the 
implications of our analysis for the viability of the CDM 
scenario, the discrepancy between X-ray and lensing masses 
in Abell 1689, and the presence of non-thermal pressure. 

3.1 Best-fit parameters 

In Figure [T] we present the results of our joint analysis for 
Abell 1689. Both the X-ray and lensing data are reasonably 
fitted by our model, with a total Xtot = Xx + Xsi + Xwi = 33.8 
(with 20 degrees of freedom). The separate contributions are 
X 2 = 12.4, 6.3, 15.1 for X-ray, SL and WL data, respectively, 
where xLs « -21og(£i en8 )- 

Note that in Figure [T] the WL model predictions are 
slightly in disagreement with the data, overpredicting the 
surface mass profile in the range 500-2000 kpc at around the 
1-cr significance level, while the X-ray and SL data are rather 
well fitted. We checked that the inclusion of WL data does 
not significantly change the values of the desired physical pa- 
rameters. Also, we verified that the joint fit is preferentially 
more strongly constrained by the X-ray and SL data with the 
present datasets, WL measurements being of poorer quality. 
Excluding WL data slightly increases (5%-15%) the errors 
on the final parameters with respect to the case where we in- 
clude them, with the exception of P nt , where the error would 
increase of ~ 60%. With this respect, while for the present 
datasets WL data are less relevant in constraining the model 
parameters with respect to X-ray and SL data, whose qual- 
ity is exceptionally good, we observe that this might not 
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Figure 1. Joint X-ray and lensing analysis for Abell 1689. In the upper panel we display the two quantities that enter in the X-ray 
analysis (Equation l llOl l): the observed projected temperature T* ro ^ . (big points with error bars) and the theoretical projected temperature 
7 1 proj,i(<j) (diamonds). We also show the theoretical three-dimensional temperature Tj(q) (points), which generates T pro j i(q) through 
convenient projection techniques. In the middle and lower panel we display the two quantities which enter in the WL and SL analysis, 
respectively (Equation lilt ): the observed surface mass profile EJ (points with error bars) and the theoretical one S(q) (solid line). The 
distance of the points for both Tj (q) and S* is representative of the true spatial resolution of X-ray and lensing data, respectively. 



apply for other galaxy clusters and other datasets, where 
the inclusion of WL data could be important in removing 
the degeneracy among the same parameters. 

To further probe the disagreement between WL and 
SL data, we jointly fitted them (without the X-ray data) 
with a spherical NFW (a = 1) model first (following 
IfJmetsu fc Broadhurstl |2008)), and then with the more 
general triaxial gNFW. We found a marginal disagree- 
ment between SL and WL data for both models. A sim- 
ilar slight disagreement can also be seen in Figure 13 of 
lUmetsu fc Broadhurstl |2008T l. 

In order to make a comparison with Mor andi et al.l 
|201ll ). we also compared the \ 2 by neglecting WL data 
and wit hout P n t ■ we fin d \ 2 = 25.3 (14 degrees of freedom), 
while in iMorandi et alj i|201ll ) we found that the \ 2 ~ 7.4 
(11 degrees of freedom). Note that the difference in the value 
of the x 2 reflects a different treatment of the SL data in the 
present paper ( i)2.2l and i)2.3p in comparison to our previ- 
ous work. Among these difference, the major improvement 
is based on the determination of covariance C s \ among the 
pixels of the 2D map distribution, in order to infer the co- 



variance C s \ among the rebinned measurements (Equation 
[3]l , and on an improved rebinning schema (Equation [2]| . We 
also point out that this difference in the value of the \ 2 
translates to slightly different errors in the parameter val- 
ues, with negligible impact on their expectation values. For 
example, we find that the axis ratio along the line of sight 
Vdm.c' = 2.39 ± 0.02 (by neglecting WL data and without 
P n t), whereas the value deduced in IMorandi et all (|201lD 
was 77dm, c' = 2.37 ± 0.11. 

In Table [T] we present the best-fit model parameters for 
our analysis of Abell 1689. Errors in the individual parame- 
ters (q, o- S ys) have been evaluated by considering the average 
value and the mean absolute deviation of the marginal prob- 
ability distribution of each parameter. Note that we list the 
axis ratio along the line of sight ?7rjM,o' rather than the eccen- 
tricity. Our work shows that Abell 1689 is a triaxial galaxy 
cluster with DM halo axis ratios ?7DM,b' = 1-24 ± 0.13 and 
?7dm.c' = 2.02 ± 0.01, where r? DM ^ b / is the axis ratio on the 
plane of the sky inferred from SL measurements, and r?DM,c' 
is the axis ratio along the line of sight inferred through our 
joint analysis. Note that these elongations are statistically 
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Figure 2. Joint probability distribution of r/ D M c ' an d £■ The 
solid(dashed) line represents the 1(2)-<t error region, while the 
big dot represents the best-fit values. 



significant, i.e., our results disprove the spherical geometry 
assumption. The axis ratio of the gas is r? ga s,b' ~ 1.1 — 1-06 
(on the plane of the sky) and r/ g ^ s c i ~ 1.5 — 1.3 (along 
the line of sight), moving from the center toward the X-ray 
boundary. 

In Table [T] we also report the value of M200 : 



8007T , . . , . , 3 

M200 = — - — (a r /c r ){br/c r )H2oo p c 



(13) 



where p c is the critical density of the Universe at the redshift 
2 of the cluster, c = i?20o/?*s, and a r /c r and b r /c r are the 
minor-major and intermediate-major axis ratios of the DM 
halo, respectively. 

The second main result from our work is the need for 
a appreciable non-thermal pressure support, formally at a 
level of ~20%. 

In Figure [2] we present the joint probability distri- 
bution of 77dm,c' an d £■ Note the positive correlation be- 
tween £ and 77dm, c'i i- e - the X-ray/Lensing mass discrep- 
ancy in clusters with prominent strong lensing features can 
be explained via a combination of both triaxiality and non- 
thermal support. This expected in light of the following con- 
siderations: (1) the observed temperature profile and pro- 
jected mass profile are both sensitive to triaxiality, specifi- 
cally T(R, r? DM ) ~ tw c ' T(R, i/dm = 1) and E(R) oc f? DM ,c' 
l|Morandi et al J [20101 ). so that the dependency of E(i?) on 
triaxiality is stronger than that of T(R) (remembering that 
generally ?7dm,c' > Vgs. B ,c')', (2) the non-thermal pressure 
support affects only the X-ray data via Equation 21 i.e., 
neglecting P nt systematically lowers the determination of 
cluster masses based on X-ray data roughly by the factor £. 



3.2 Implications for the viability of the CDM 
scenario and the X-ray/Lensing mass 
discrepancy 

A precise determination of cluster masses is crucial for the 
use of clusters as cosmological probes. However, there are 
discrepancies between cluster masses determined based on 
gravitational lensing and on X-ray observations, the former 
being significantly higher than the latt er in many clusters 
with p rominent lensing features. Indeed. [Oeuri fc Blandfordl 

(2009) showed that SL clusters with the largest Einstein 
radii constitute a highly biased population with major axes 
preferentially aligned with the line of sight thus increas- 
ing the magnitude of lensing. Given that lensing depends 
on the integrated mass along the line of sight, either for- 
tuitous alignments with mass concentrations that are not 
physically related to the galaxy cluster or departures of 
the DM halo from spherical symmetry can bias upwards 
the three-dime nsional mass with respect to the X-ray mass 
|Gavazzill20"05l ); on the other hand, X-ray-only masses hinge 
on the accuracy of the assumption of strict hydrostatic equi- 
librium: the presence of bulk motions in the gas can bias 
low the three-dimension al mass profile between 5 and 20% 
jMeneghetti et alJl2oToT ). 

lOeuri et all l|2005t ) concluded that weak lensing mea- 
surements in Abell 1689 are indeed compatible with the 
CDM-based triaxial halo model if Abell 1689 represents a 
rare population (~6% b y number) of cluster-scale halos, 
and lMorandi et all (|201ll ) demonstrated that triaxiality al- 
lows us to remove the mass discrepancy betw een the strong 
lensin g and X-ray estimates in Abell 1689. iMolnar et al.l 

(2010) suggested an alternative explanation for the mass 
discrepancy, due to violations of strict hydrostatic equilib- 
rium in Abell 1689: they found that a contribution of about 
40% from non-thermal pressure within the core region of 
Abell 1689 can explain the mass discrepancy, provided that 
the spherical geometry assumption holds. 

Moreover, recent work investigating the mass distribu- 
tions of individual galaxy clusters (Abell 1689 and oth- 
ers) based on gravitational lensing and employing stan- 
dard spherical modeling have found a potential inconsis- 
tency compared to the predictions of the CDM scenario re- 
lating halo mass to the concentration parameter c. In partic- 
ular, relatively high values of c (~ 8-14) have been derived 
from lensing analysis of Abell 1689 (|Broadhurst et aTll2005l : 
iLimousin et al.l feoO'it ). These values are outside the range 
predict ed from simulat ions of the standard CDM model 
(c ~ 4; iNeto et alj|2007h . 

In the present work we presented a physical model for 
the cluster Abell 1689 with a triaxial mass distribution in- 
cluding support from non-thermal pressure. This model is 
consistent with X-ray and SL observations as well as with the 
predictions of the CDM model. Doing so we have removed 
the apparent discrepancy between X-ray and lensing mass 
estimates, and potential inconsistencies discussed in the lit- 
erature between the predictions of the CDM scenario and 
the measurements in clusters with prominent strong lensing 
features. One of the main results of our work is the mea- 
surement of a central slope of the DM profile 0.95 ± 0.05 in 
agree ment with the theore t ical expectation of the CDM sce- 
nario (INavarro et al. Hl997l : lMerriU et alj|2006l ). and a value 
of the concentration parameter 5.71 ± 0.47, in agreement 
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Table 1. Best-fit model parameters of A1689. The columns 1—6 refer to the best fit parameters c,r s ,a,r] DMtC i , £ and <T sys , while the 
last column refers to M200- 



c r s 


a VDM,c> 






M 20 o 


(kpc) 






(g cm -2 ) 


(1O 15 M ) 


5.71 ± 0.47 348 ± 28 


0.95 ±0.05 2.02 ±0.01 


0.22 ± 0.01 


0.008 ± 0.001 


2.59 ±0.09 



wi th the theo r etical expectation from N-body simulations 
of INeto et all l|2007l ), where c ~ 4 at the redshift and for 
the virial mass of Abell 1689, with an intrinsic scatter of 
~ 20%. In this way we also managed to reproduce the large 
Einstein radius of A bell 1689 via our triaxial framework 
(|Morandi et al.ll201ll ). 

Even with triaxiality, though, lOeuri et al l {2005) found 
that Abell 1689 is relatively rare in its lensing efficiency , 
as noted above. Furthermore, Broadhurst & Bar kanal (|2008l ) 
showed that the observed Einstein radius #e of four well- 
studied massive clusters (including Abell 1689) lies well be- 
yond the predicted distribution of Einstein radii in the stan- 
dard CDM model, typically by a factor of 2, even after ac- 
counting for the actual projected mass distributions in N- 
body simulations, and for lensing bias. Using a large r sam- 
ple of 12 X-ray selected clusters, IZitrin et al.l (l201ll ) found 
a smaller discrepancy with CDM predictions, typically by a 
factor of 1.4. 

In a similar fashion, by comparing numerical and obser- 
vational cluster samples, [Mencghc tti et ail (|201ll ) show that 
some real clusters have too large lensing cross sections and 
Einstein rings compared to expect ations in a AC DM cosmo- 
logical model. iHoresh et alj (|201ll ) partially confirmed these 
results, by comparing the lensed arc statistics measured from 
the Millennium simulation to those of a sample of observed 
clusters: they found an excellent agreement between the ob- 
served and simulated number of arcs in the redshift range 
0.3 < z < 0.5, while at lower redshift some conflict still re- 
mains, with real clusters being ~ 3 times more efficient arc 
producers than their simulated counterparts. This departure 
might be due to selection biases in the observed subsample 
at this redshift or to physical effects that arise at low red- 
shift and enhance the lensing efficiency, and not included 
the simulat ions, as the effect of b aryons on the DM profile 
in clusters l|Barkana fc Loebll2010l ). 

We also compare the major-minor principal axis ratio 
ijom^ = 2.02 ±0.01 wi th that inferred in ou r previous work 
(Vdm,c> = 2.37 ±0.11; iMorandi et~al"1l201ll ). The estimated 
value of the elongation along the line of sight in the present 
work is somewhat smaller, due to the assumption of strict 
hydrostatic equilibrium and the lack of the large-scale WL 
data in our previous work. This lends further support to our 
emphasis on the role of the effects of both geometry and non- 
thermal pressure support on the physical parameters and 
open a new window in recovering the intrinsic shapes and 
desired physical parameters of galaxy clusters in a bias-free 
way. The estimated value of 77dm, c' is also consistent w ith 
the results from numerical simulations (|Shaw et al.l l2006). 



3.3 Non-thermal gas pressure 

Simulations of galaxy clusters predict that turbulent mo- 
tions should occur in the IC gas while the matter contin- 
ues to accrete along filaments. This energy should then cas- 
cade from large to small scales and can eventually dissipate 
into the gas. The attainment of HE does not by itself rule 
out an appreciable level of turbulence (~5-15%, lLau et ail 
2009; [Zhang et al.ll2010l ). Measurements of the non-thermal 
energy in the IC gas are important in order to estimate 
the amount of energy injected into clusters from mergers, 
accretion of material or feedback from active galactic nu- 
clei (AGNs); in the latter case, this probes energy trans- 
port from th e central nucleu s into its surroundings and into 
the IC gas. iBriiggen et "all (|2005l ) found that in clusters 
with AGN feedback gas motion induced by the inflation of 
bubbles and their buoyant rise leads to velocities of about 
500-1000 km s" 1 , which is a significant fraction of the local 
sound speed. Strong turbulent motions would provide sig- 
nificant non-thermal pressure support; this would bias low 
the determination of cluster mass profiles measured under 
the assumption that P to t = Pth in Equation Q. The mea- 
surement of turbulence is thus important for mass determi- 
nation as wel l as un derstanding AGN feedback. Moreover, 
IMolnar et aTJ (|2010| ) argued that there is no need for feed- 
back from a central AGN for the strict hydrostatic equilib- 
rium to break down in the central regions of clusters: sub- 
sonic random gas motions, a direct consequence of hierar- 
chical structure formation, allow us to explain violations of 
the strict hydrostatic equilibrium assumption in the central 
regions; in the outer regions violations could be due to more 
recent and still ongoing slow accretion, since those regions 
have not reached equilibrium due to the large sound crossing 
time. 

Observationally, IRichard et all (|2010l ) measured an X- 
ray/Lensing mass discrepancy of < Af s i/Af x >= 1.3 at 3- 
cr significance level in a sample of 20 strong lensing clus- 
ters. They interpret this as evidence that the assumption of 
strict hydrostatic equilibrium required by the X-ray mass es- 
timates is not wholly reliable and the merging activity can 
add non-ther mal pressure support to the IC gas through 
bulk motions. ISanders et alj ([2010) placed a direct limit on 
turbulence based on the non-thermal velocity broadening 
measured from the emission lines originating in the central 
30 kpc of the galaxy cluster Abell 1835. They found that 
the ratio of turbulent to thermal energy density in the core 
is less than 13%. 

IMolnar et al.l (|2010l ) analyzed a sample of massive clus- 
ters of galaxies drawn from high-resolution cosmological sim- 
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ulations and found a significant contribution (20%-45%) 
from non-thermal pressure. They also tested the valid- 
ity of strict hydrostatic equilibrium in Abell 1689 using 
gravitational lensing and X-ray observations under the as- 
sumption of spherical geometry in order to explain the X- 
ray /Lensing mass discrepancy: they found a contribution of 
about 40% from non-thermal pressure within the core re- 
gion of Abell 1689, suggesting an alternate explanation for 
the mass discrepancy, as long as the spherical assumption 
holds. 

Our results point to a scenario where the non-thermal 
component is about 20% of the total energy budget of 
the IC gas. This lev el is quite lower than that found from 
iMolnar et ail (|2010h under the assumption of spherical ge- 
ometry, suggesting that accounting for the proper triaxial 
geometry is quite important in evaluating P nt accurately. 

Another relevant consideration here is that in our model 
we constrained the non-thermal pressure to be a constant 
factor of the loca l thermal pressure th roughout the clus- 
ter. The results of iMahdavi et all (2008) show that there is 
a radial trend of the X-ray/WL mass ratio, that is inter- 
preted as caused by non-thermality increasing toward the 
outer regions, though their findings hinge on the assumed 
spherical geometry, so they did not disentangle the effect of 
triaxiality from non-thermal pressure support. In this per- 
spective, we used a linear relation of P nt , by fixing the slope 
to that found bv lMahdavi et all l|2008h and leaving the nor- 
malization as a free parameter: the results are not affected 
appreciably from this, the reason is that the fit is more domi- 
nated by the innermost part (< 300 kpc) of the Xray and SL 
data. Moreover, in order to gauge the likely model depen- 
dence of the inferred non-thermal pressure, we performed 
a joint analysis by excluding the temperature constraints 
(Equation [TUJ) and assuming a spherical NFW model for the 
DM. Doing so also allow s us to make a direct comparison 
with lLemze et all ((2008), who used a similar approach. We 
found that the ratio between the (de)-projected tempera- 
ture profile measured from the X-ray spectrum and that 
obtained by assuming strict hydrostatic equilibrium (Equa- 
tion U with Ptot = Pth) is ~ 0.5 — 0.6, which is lower than 
the deter mination in our tri axial model (Pth/Aot ~ 0.8). 
Note that|L emze et al.l (|2008l ) concluded that a factor of 0.7 
explains most of the temperature discrepancy in Abell 1689, 
whereas we find a somewhat lower value; this stems mostly 
from the more recent Chandra calibrations adopted in the 
present work, which led to a lower spectral temperature pro- 
file (by ~ 20%). While this test heavily hinges on the model 
of spherical geometry, it suggests that the small formal error 
on £ in our triaxial joint analysis merits some caution. 



4 SUMMARY AND CONCLUSIONS 

In this paper we have employed a physical cluster model 
for Abell 1689 with a triaxial mass distribution including 
support from non-thermal pressure, and proved that it is 
consistent with the X-ray and SL observations and the pre- 
dictions of CDM models. 

We demonstrated that accounting for the three- 
dimensional geometry and the non-thermal component of 
the gas allows us to resolve the long-standing discrepancy 
between the X-ray and strong lensing mass of Abell 1689 in 



the literature, as well as to measure a central slope of the 
DM and a concentration parameter in agreement with the 
theoretical expectations of the CDM scenario. 

We also measured the contribution of the non-thermal 
component of the gas (~ 20% of the total energy budget of 
the IC gas). This has important consequences for estimating 
the amount of energy injected into clusters from mergers, 
accretion of material or feedback from AGN. 

The increasing precision of observations now makes it 
possible to test the assumptions of spherical symmetry and 
hydrostatic equilibrium. Since important current cosmolog- 
ical tests are based on the knowledge of the masses, shapes, 
and profiles of galaxy clusters, it is important to better char- 
acterize their physical properties by allowing for realistic 
triaxial structures as well as non-thermal pressure support. 
The application of our method to a larger sample of clusters 
will allow to infer the desired physical parameters of galaxy 
clusters in a bias-free way, with important implications on 
the use of galaxy clusters as precise cosmological probes. 
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